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MOLECULAR  BEAM  CHEMISTRY:  REACTIONS  OF  OXYGEN  ATOMS  WITH  HALOGEN  MOLECULES 

+  * 

D.  D.  Parrish,  D.  A.  Dixon,  and  D.  R.  Herschbach* 

Department  of  Chemistry,  Harvard  University 
Cambridge,  Massachusetts  02138 

Molecular  beam  experiments  show  that  the  reaction  0  +  -*•  BrO  +  Br 

proceeds  via  formation  and  decomposition  of  a  collision  complex  which  persists 

for  many  vibrational  periods  and  at  least  a  few  rotational  periods,  or 
'v  -12 

>5  x  10  sec.  The  angle  and  translational  energy  distributions  agree  closely 
with  those  predicted  by  approximate  statistical  complex  models  akin  to  the  RRKM 
theory  of  unimolecular  decay.  The  analogous  reactions  of  I-,  IC1  and  I3r 
were  also  found  to  proceed  via  long-lived  complexes.  For  IC1  and  IBr  a  large 
yield  of  10  was  found,  but  no  CIO  or  BrO,  despite  a  more  favorable  reaction 
exoergicity  for  these  latter  products.  These  results  suggest  the  intermediate 
complexes  are  unsymmetrical  species,  0-Br-Br  or  0-I-C1,  which  are  probably 
stable  by  >20  kcal/mole  and  in  triplet  spin  states.  The  reaction  0.  +  C^  is 
found  to  proceed  most  likely  via  an  osculating  complex.  The  reactive  surface 
does  not  sample  the  well  corresponding  to  C2v  0%2  in  its  ground  singlet  state. 
The  complex  is  predicted  to  be  a  linear  or  slightly  bent,  asymmetric  triplet 
0XX  with  an  electron  configuration  of  either  •••(o*)  (tt  )  or  •••(tt  )  (a  )  . 


^Present  address:  Metropolitan  State  College,  Denver,  Colorado 

^Present  address:  Chemistry  Jpepartment,  University  of  Minnesota,  Minneapolis,  MU. 


INTRODUCTION 


The  electronic  structure  of  reactant  and  product  molecules  and  transition 
states  governs  reaction  kinetics.  While  much  has  been  inferred  from  macrosco¬ 
pic  experiments,  only  single  collision  studies  can  hope  to  examine  the  actual 
molecular  dynamics.*  Even  such  studies  may  not  be  able  to  abstract  detailed 
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information  about  the  potential  surface.  However,  significant  qualitative 
electronic  structure  information  about  transition  states  or  reaction  interme¬ 
diates  can  often  be  obtained.  In  favorable  cases,  it  is  possible  to  synthesize 
intermediates.  For  reactions  that  proceed  through  a  long-lived  complex,  even 
qualitative  dynamics  can  yield  information  about  the  stability  of  these  systems 
that  is  not  presently  available  from  bulk  experiments  or  from  electronic 
structure  calculations. 

This  paper  reports  a  study^  of  several  reactions  of  oxygen  atoms  with 
halogen  molecules, 

0(3Pg)  +  XY(*Z+)  -  0X(2tt)  +  X(2PU),  (Rl) 

with  XY  =  0*2 »  Br2,  I2>  I CJt ,  and  IBr.  Some  of  these  reactions  have  been  ex- 
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amined  by  classical  kinetic  methods  using  flow  tubes  and  diffusion  flames, 
and  concurrent  molecular  beam  experiments  have  been  pursued  in  other  labora- 
tories.  *  We  find  that  at  collision  energies  of  2-3  kcal/mol  these  reactions 
proceed  via  a  long-lived  collision  complex,  except  for  the  Ci^  case  which 
proceeds  via  a  osculating  complex.  At  higher  collision  energies  there  is  now 
evidence  that  the  other  0  +  XY  complexes  likewise  become  osculating. ^  Since 
the  ratio  of  the  reduced  mass  of  the  reactants  to  that  of  the  products  is 
rather  small,  the  0  +  XY  reactive  scattering  proves  to  be  particularly  suita¬ 
ble  for  comparison  with  simple  statistical  models.  These  treat  the  complex 
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mechanism  by  phase  space  and  transition  state  models  patterned  after  the 
RRKM  theory  of  unimolecular  decomposition.** 


The  analogous  reactions  of  hydrogen  atoms,  methyl  radicals,  alkali  atoms 


alkali  dimers,  alkaline  earth  atoms,  halogen  atoms,  and  mercury  atoms  with 

12 

halogen  molecules  have  also  been  studied  in  molecular  beams.  This  family 

of  reactions  scans  the  full  panoply  of  reaction  dynamics,  from  extremely  im- 
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pulsive  processes  akin  to  photodissociation  in  the  H  +  XY  case  to  statis¬ 
tical  behavior  akin  to  unimolecular  decomposition  via  a  well  defined  interme- 
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diates  in  the  Hg  +  XY  case.  The  0  atom  reactions  have  a  special  place  in 
this  family  because  only  they  offer  the  possibility  of  transitions  between 
triplet  and  singlet  potential  energy  surfaces.  In  (Rl),  the  reactants  approach 
on  a  triplet  surface,  whereas  the  products  can  depart  on  either  a  singlet  or 
triplet  surface.  In  accord  with  qualitative  electronic  structure  arguments,4 
the  reactive  scattering  indicates  that  these  reactions  go  predominantly  via  a 
triplet  0-X-Y  surface  without  transition  to  the  singlet  X-O-Y  surface  corres¬ 
ponding  to  the  ground  electronics  state  of  the  known  halogen  oxide  molecules. 
This  result,  interpreted  in  terms  of  orbital  correlations  and  electronegativity 
differences,  led  to  the  prediction4  that  the  0  +  reaction  would  require  a 
relatively  large  activation  energy  even  though  more  exoergic  than  the  other 
0  +  XY  reactions.  This  was  subsequently  confirmed  in  a  flow  tube  study  ,  which 
found  Eact  =  12  kcal/mol  for  the  0  +  case,  in  contrast  to  EflC^  =  2  and  zero 
kcal/mol  for  the  Cl 2  and  Br~  systems. 


EXPERIMENTAL 

mvvvvwvw/ 


The  apparatus  is  schematically  depicted  in  Fig.  1.  Only  key  features  and 
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recent  modifications  will  be  described  here;  more  details  are  given  elsewhere  ’  . 

The  two  reactant  beams  cross  at  an  angle  of  90°  in  a  large  scattering  chamber.  A 
4200  1/sec  oil  diffusion  pump  and  a  large  liquid  nitrogen  cooled  shield  pump  non¬ 
condensible  and  condensible  gases,  respectively,  in  this  chamber.  Reaction  products 
formed  in  the  col  1 ision  zone  are  detected  by  a  mass  spectrometer  mounted  on  the 
lid  of  the  scattering  chamber.  This  lid  is  rotatable  allowing  the  detector  to 


be  moved  In  the  plane  of  the  beams  to  scan  the  angular  distribution  of  the 
scattered  reaction  products.  The  detector  Is  comprised  of  an  electron-bombard¬ 
ment  Ionizer,  quadrupole  massfilter,  scintillation  Ion  counter,  and  two  gated 
scalars  synchronized  with  the  modulation  of  the  halogen  molecule  beam.  In 
the  velocity  analysis  experiments,  a  slotted  disk  chopper  is  attached  to  the  lid 
between  the  collision  zone  and  the  detector.  A  PDP  8/L  computer  is  used  both 
as  a  multiscalar  to  record  the  time  of  flight  (TOF)  spectrum  of  the  products 
and  for  immediate,  approximate  extraction  of  the  velocity  distribution  from 
the  recorded  spectrum. 

To  reduce  the  residual  partial  pressure  of  background  product  molecules 
in  the  ionization  region  of  the  detector  to  a  low  level  (£l(f^  torr)  the 
ionizer  is  nested  in  three  chambers  differentially  pumped  by  50  1/sec  ion  pumps. 
The  ionizer  is  surrounded  by  liquid  nitrogen  cooled  surfaces  and  the  design 
allows  unionized  product  molecules  (^99.9%  of  the  total)  to  "fly  through"  the 
ionizer.  A  new  pump  element  which  is  a  46  cm  x  2.5  cm  x  0.32  cm  gold  plated 
copper  strip  cooled  to^l^K  by  an  Air  Products  and  Chemicals  Inc.,  Displex 
Model  CS-202  closed  cycle,  displacer/expanded  helium  refrigerator  has  been 
added.  The  l^K  element  is  completely  surrounded  by'W/^K  radiation  shields 
cooled  by  the  refrigerator' s  first  stage.  The  refrigerator  is  mounted  at  the 
throat  of  the  ionizer  chamber  ion  pump  and  the  pumping  element  extends  down  this 
chamber  to  within  13  cm  of  the  ionizer.  This  refrigerator  has  given  more  than 
6000  hours  of  maintenance  free  operation. 

Beam  Sources.  The  oxygen  atom  beam  is  produced  by  a  low  pressure  (£l00 

microns  Matheson  Extra  Dry  oxygen)  microwave  discharge.  The  source  shown  in 
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Fig.  2a  is  similar  to  that  used  previously  for  deuterium  atoms  and  chlorine 
atoms^®.  The  discharge  cavity  is  essentially  the  foreshortened  1/4  wave  coaxial 
cavity  described  by  Broida  and  coworkers  .  However,  for  convenience  of  mounting 
and  vacuum  connections,  a  wave  guide  has  been  added  by  lengthening  the  inner  and 
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outer  electrodes  of  the  cavity  by  four  wavelengths.  The  discharge  occurs  In  a 
quartz  tube  which  Is  sealed  Into  the  gold  plated  copper  waveguide  with  epoxy 
resin  to  provide  a  vacuum  tight  seal.  Water  cooling  colls  protect  the  epoxy 
from  pyrolysis  while  air  cooling  prevents  melting  the  quartz  tube.  Two 
important  requirements  for  satisfactory  operation  have  been  identified:  the 
quartz  tube  must  be  coated  with  freshly  prepared,  syrupy  phosphoric  acid  before 
each  run,  and  the  dimension  indicated  by  A  in  Fig.  2a  must  be  held  as  small  as 
possible.  The  microwave  frequency  is  2450  MHz  and  is  operated  at  'v/lOO  watts 
power.  The  beam  exit  slit  (1  mm  x  10  mm)  is  located  ^2  cm  from  the  center  of 
the  cavity.  Since  the  beam  is  taken  from  the  region  of  active  discharge  it 
could  be  expected  to  contain  significant  amounts  of  0 ^  (U)  and  0(^D)  metastable 
species  in  addition  to  ground  state  0(  P)  atoms.  The  A)  is  unreactive  at 
the  experimental  collision  energy  with  the  halogens  as  all  reactions  are  at 
least  16  kcal/mole  endoergic.  The  0(^0)  is,  however,  generally  quite  reactive 
and  would,  if  present,  possibly  interfere  with  the  observation  of  the  products 
of  ground  state  0  atom  reactions.  However,  the  state  is  efficiently  collision- 
ally  quenched  by  0^  and  is  thus  not  expected  to  be  a  significant  component  of 
the  beam;  indeed,  we  have  seen  no  evidence  for  reactions  of  this  excited  species 
as  is  discussed  later.  Velocity  analysis  of  the  beam  indicate  that  it  is 
trans  latior.ally  hot  (1000°  to  1200°K)  with  a  full  Maxwellian  velocity  distribution. 
The  fluxes  of  0  and  0^  are  approximately  equal  under  our  optimum  discharge 
conditions. 

The  halogen  molecule  source  shown  in  Fig.  2b  is  a  high  pressure  (150  to 
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900  torr)  supersonic  nozzle  .  The  nozzle  is  fabricated  from  a  6.2  mm  0.D  quartz 

// 

tube  pulled  down  to  form  a  short  closed  tip.  This  tip  is  ground  back  until  a  0.10 
to  0.17  mm  dia  orifice  is  produced.  A  heating  mantle  around  the  quartz  tube  allows 
tne  nozzle  to  be  operated  from  300°  to  1000°K.  The  nozzle  is  moveable  along  the 
beam  axis  and  is  mounted  at  an  optimum  value  of  mm  behind  a  1.5  mm  skimmer. 


The  nozzle  beam  chamber  is  pumped  by  a  large  liquid  nitrogen  cooled  shield  and 
a  2400  1/sec  oil  diffusion  pump  to  achieve  the  very  high  pumping  speeds  required. 
The  resolution  of  the  velocity  analyzer  is  insufficient  to  verify  the  very 
narrow  velocity  distributions  expected  for  this  beam.  The  velocity  analysis  is 
consistent  with  high  Mach  numbers  {>10)  and  nearly  complete  relaxation  of  the 
rotational  degrees  of  freedom  of  the  halogen  molecules. 

Both  beam  sources  are  mounted  in  chambers  differentially  pumped  by  1200  1/sec 
oil  diffusion  pumps.  Collimating  slits  (1.25  mm  x  10  mm)  lead  to  the  scattering 
chamber.  The  measured  angular  widths  of  the  resulting  beams  are  5°  (full  width 
at  half- peak  intensity).  Detailed  calculations  indicate  this  is  the  maximum 
size  consistent  with  a  scattering  region  cross-section  small  enough  to  be 
completely  viewed  by  the  detector.  A  rotating  paddle  modulates  the  halogen  beam 
at  70  Hz  and  a  beam  flag  is  used  to  interrupt  the  oxygen  atom  beam.  The  beam 
modulator  and  the  beam  flag  are  located  in  the  respective  differentially  pumped 
chambers. 

Experimental  Procedures.  The  angular  distributions  are  obtained  (with  the 
TOF  chopper  removed)  by  counting  at  each  laboratory  scattering  angle'  for  a  period 
of  20  to  200  sec  with  the  beam  flag  open  and  then  iimiediately  for  the  same  period 
with  the  beam  flag  on  the  0  atom  source  closed.  The  modulated  signal  (difference 
between  "on”  and  "off"  scalars)  with  the  beam  flag  closed  is  subtracted  from  the 
open  beam  flag  modulated  signal  to  obtain  the  reactive  scattering  intensity.  In 
the  studies  reported  here  the  closed  beam  flag  signals  were  typically  0  to  30% 
of  the  beam  flag  open  signals  and  were  angular  dependent,  increasing  to  a  maximum 
near  the  chopped  halogen  beam;  these  beam  flag  closed  signals  can  be  attributed  to 
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the  reaction  of  the  halogen  molecules  with  background  oxygen  atoms  in  the  scatter¬ 
ing  chamber.  The  backgrounds  ("off"  scalar)  were  1  to  20  times  larger  than  the 
reactive  signals.  To  obtain  the  angular  distributions  four  separate  scans  through 
the  complete  angular  range  were  made  in  varying  angle  sequence.  Since  each  data 
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point  was  relatively  noisy,  the  time  variation  in  the  signal  at  each  angle  was 
included  in  the  procedure  developed  to  time  normalize  for  beam  Intensity  or 
detector  sensitivity  variations  during  the  experiment.  This  method  is  described 
in  the  Appendix.  It  gives  significantly  better  data  than  the  traditional  method 
of  using  the  time  variation  of  the  signal  at  a  single,  relatively  noise-free 
angle  for  normal ization.  The  time  normalized  signals  were  averaged  over  the  four 
angle  scans  to  obtain  the  data  points.  Simple  error  propagation  techniques  are 
described  in  the  Appendix  which  allowed  the  systematic  derivation  of  confidence 
limits  and  error  bars. 

The  TOF  velocity  analysis  data  are  collected  in  the  manner  described  pre¬ 
viously.  A  slightly  narrower  chopper  disk  slit  width  (0.15  cm)  yielded  a 
shutter  function  of  0.067  rc  (full  width  at  half  maximum)  with  a  chopper  period 
Tr  of  860  usee  and  channel  width  of  24  usee.  Some  improvements  have  been  made  in 
the  procedure  for  converting  the  TOF  spectrum  to  the  velocity  distribution.  At 
low  velocities  (long  flight  times),  the  velocity  resolution  (Av/v^  with  Av 
the  FWHM  at  nominal  velocity  v)  is  determined  primarily  by  the  ratio  of  ionizer 
length  (2.5  cm)  to  length  of  flight  path  (13.5  cm)  and  is  relatively  unaffected 
by  the  TOF  channel  width.  Hence,  two  or  more  consecutive  channels  at  the  end  of 
the  TOF  spectrum  can  be  combined  by  adding  the  counts  in  each;  the  analysis  with 
the  resulting  longer  channel  width  gives  improved  signal-to-noise  ratios  with 
only  a  small  decay  in  velocity  resolution.  In  practice,  the  data  were  analyzed 
by  beginning  with  the  last  channel  and  combining  with  preceding  channels  sequen¬ 
tially  until  the  product  of  the  velocity  resolution  and  the  noise-to-signal  ratio 
was  minimized.  This  process  weys  then  repeated  with  each  channel.  Error  propagation 
techniques  again  allowed  assignment  of  systematic  error  bars.  The  low  velocity 
data  exhibit  less  scatter  than  would  be  expected  from  the  magnitude  of  the  error 
bars.  This  phenomenon  arises  because  consecutive  points  come  from  overlapping 


time  periods  and  hence  Include  much  of  the  same  data.  The  error  bars  do,  however. 
Indicate  the  random  errors  in  the  data. 

RESULTS 

Table  I  gives  the  important  parameters  that  varied  between  experimental  runs. 
For  a  given  reactant  system,  identical  conditions  were  maintained  as  closely  as 
possible  for  the  angular  velocity  distribution  runs.  The  measured  product  angular 
distributions  are  shown  in  Figs.  3  and  4  along  with  curves  calculated  from  the 
statistical  complex  models  discussed  later.  In  all  cases  the  angular  distribution 
peaks  near  the  direction  of  the  nominal  centroid  vector. 

The  0  +  Cl g  distribution,  however,  does 
exhibit  a  slight  shift  toward  the  0  atom  beam.  These  experimental  angular 
distributions  exhibit  no  clear  indication  of  the  bimodal  structure 

generally  associated  with  a  collision  complex  mechanism.  This  lack  of 
structure  is  due  to  the  inability  of  the  experiments  to  resolve  the  two  peaks 
which,  for  these  systems,  are  very  close  together  in  angular  range.  As  discussed 
below  this  behavior  is  clearly  understood  and,  indeed,  is  directly  predicted  by 
the  collision  complex  model. 

Careful  searches  were  made  for  the  other  possible  product  (CIO  from  0  +  IC1 
and  BrO  from  0  +  IBr)  from  the  mixed  halogen  reactions,  but  no  indications  of  these 
products  were  found.  Energetic  considerations  (Table  II)  do  not  prohibit  the 
formation  of  these  alternative  products  and  statistical  phase  space  theory  predicts 
that  CIO  should  be  the  dominant  product  of  the  0  +  IC1  reaction.  Figure  5  indicates 
the  sensitivity  of  the  detector  for  the  CIO  product  and  shows  that  the  signal  from 
this  product  must  be  at  least  one^  order  of  magnitude  less  intense  than  expected 
from  the  statistical  complex  model  assuming  a  distribution  of  products  from  phase 
space  theory.  The  low  intensity  of  the  signal  from  the  0  +  IBr  reaction  prevented 
a  definitive  search  for  the  BrO  product  from  this  reaction. 
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The  measured  product  velocity  spectra  are  shown  in  Figs.  6-8  along  with 
curves  calculated  from  the  statistical  complex  models  for  the  reactions  0  +  Cl2t 
0  +  Br2  and  0  +  IC1.  Signals  were  so  small  for  the  0  +  Cl ^  reaction  (the  single 
velocity  distribution  required  3  hr.)  that  only  one  velocity  distribution  near 
the  peak  of  the  angular  distribution  was  measured.  For  the  0  +  Br2  and  0  +  1C1 
reactions  each  of  the  independently  measured  velocity  distributions  are  normalized 
to  a  common  intensity  scale  by  requiring  that  the  integral  of  the  velocity 
distribution,  when  converted  to  number  densities,  equals  the  scattering  intensity 
at  the  corresponding  angle  in  the  angular  distribution  measurements. 

A  small  amount  of  background  halogen  molecules  effuses  from  the  differentially 
pumped  halogen  source  and  is,  thus,  modulated.  This  background  is  very  much  smaller 
than  the  beam  intensity  but  the  reaction  of  0  atoms  with  this  background  can 
dominate  the  detected  product  signal  at  some  ranges  of  laboratory  angles  and 
velocities.  This  effect  is  evident  at  low  laboratory  velocities  and  angles  wide 
with  respect  to  the  centroid  (6  =  102°  and  25°)  for  the  0  +  Br2  reaction.  This 
accounts  for  the  wide  deviations  from  the  model  calculations  at  these  angles; 
the  low  velocity  peak  at  6  -■  25°  could  be  attributed  to  this  reaction  with  back¬ 
ground  molecules.  Another  possible  explanation  is  due  to  the  presence  of  vjn  der 
Weals  dimers  in  the  halogen  nozzle  beam.  Some  dimers  are  expected  due  to  the 

high  source  pressures. 

The  reaction 

0  +  (Xj?  -  OXX  +  X? 

could  occur  and  the  0X2  could  fragment  in  the  ionizer  due  to  electron  bom¬ 
bardment  or  could  decay  to  OX+X  before  being  detected.  This  could  also  account 
for  the  low  velocity  peak  at  9  =  25°. 


In  all  cases  the  velocity  distributions  peak  at  velocities  that  correspond 
closely  to  forward  (0°)  or  backward  (180°)  scattering  in  the  center  of  mass 
system  and  fall  off  gradually  toward  the  energy  conservation  limits.  Poor 
velocity  resolution  of  the  products  as  well  as  energy  and  angular  spreads  in  the 
parent  beams  account  for  the  intensity  in  energetically  excluded  regions.  If 
the  parent  0  atom  beam  had  a  significant  proportion  of  metastable  0( 1 D)  atoms  which 
produced  interfering  products,  the  excess  energy  would  allow  scattering  at 
higher  laboratory  velocities.  The  lack  of  such  scattering  and  the  reasonably 
good  agreement  of  the  loose  transition  -  state  model  which  assumes  only  ground 
state  0  atoms  with  the  data  at  high  velocities  preclude  a  significant  contribution 
of  products  from  0 ( ^ D)  reactions. 

Kinematic  Analysis.  The  transformation  of  the  flux  density  I^g  in  the 
laboratory  (LAB)  system  to  the  flux  density  I™  in  the  center-of-mass  (CM)  system 
is  given  by 

ICM(e»u)  =  (uV,  ilab(0’v) 

where  u  and  0  are  the  velocity  and  angle  in  the  CM  system  and  v  and  G  the 
respective  values  in  the  LAB  system.  This  direct  transformation  with  no  averaging 
holds  true  only  for  monoenergetic  beams.  For  the  0  +  Br^  reaction  the  kinematics 
are  favorable  enough  and  the  data  are  complete  enough  to  allow  this  transformation 
to  be  made  with  reasonable  accuracy  if  the  spread  in  the  reactant  beam  velocities 
is  ignored.  This  is  accomplished  by  assigning  each  beam  its  most  probable  velocity. 
Figure  9  shows  the  results  of  this  nominal  transformation  of  the  data  for  0  +  Br? 
from  Fig.  7  with  some  smoothing  of  the  contours.  Elementary  considerations  of  the 
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collision  process  dictate  that  the  scattering  must  be  cyllndrically  symmetric 
about  the  relative  velocity  (the  0°  -  180°  axis)  and  the  data  showed  this 
within  experimental  uncertainty.  The  data  before  smoothing  showed  a  slightly 

0  Q 

greater  intensity  in  the  CM  angular  range  90  to  180  which  was  within 

22 

experimental  uncertainty.  A  reaction  which  proceeds  through  a  complex  whose 
lifetime  is  less  than  several  rotational  periods  would  be  expected  to  yield 
an  angular  distribution  with  some  asymmetry  about  0  =  90°,  but  in  such  a  case 
the  0°  to  90°  quadrant  would  have  higher  scattering  intensity,  just  opposite  to 
the  apparent  asymmetry  we  observed.  For  our  map,  we  assumed  symmetry  about 
the  0=0°  and  180°  line  and  about  0  =  90°.  These  data,  combined  with  smoothing, 
energy  conservation  considerations,  and  the  fact  that  the  contours  must  be 
perpendicular  where  they  cross  the  0°  -  180°  axis  and  the  90°  ray  give  the  map 
shown  in  Fig.  9.  The  bimodal  structure  of  the  scattering  intensity  is  clearly 
apparent  with  the  most  intense  scattering  distributed  approximately  symmetrically 
forwaro  and  backward  with  respect  to  the  initial  0  atom  direction. 

The  CM  angular  distribution  for  the  0  +  Br^  reaction  may  be  obtained  from 
Fig.  9  by  integrating  the  scattering  at  each  angle  over  the  product  velocity. 

This  result  is  presented  in  Fig.  10(a).  The  CM  velocity  distribution  at  any 
angle  may  be  obtained  by  taking  a  cut  through  the  contour  plot  of  Fig.  9  at 
the  chosen  angle.  At  all  angles  the  results  are  very  nearly  identical.  Figure 
10(b)  shows  the  CM  velocity  distribution  that  best  fits  all  angles;  Fig.  10(c) 
gives  the  product  translational  energy  distribution  that  corresponds  to  this 
velocity  distribution. 

Total  Reaction  Cross  Sections.  The  total  reactive  cross  section  could  be 
easily  derived  if  it  were  possible  to  measure  the  absolute  beam  densities  and 
the  volume  of  the  beam  intersection  region.  However,  this  determination  is  not 
possible  here  nor  has  it  been  in  any  crossed  molecular  beam  study.  Thus  indirect 
methods  are  usually  employed;  the  most  common  is  comparison  of  measured  small 
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angle  elastic  scattering  with  theoretical  predictions  to  normalize  the 

23 

experimental  Intensity  scale  to  absolute  values.  For  these  studies,  however, 
such  a  comparison  is  not  possible.  Elastic  scattering  of  0  atoms  cannot  be 
separated  from  the  mass  16  detector  background  and  mass  16  contributions  from 
the  elastic  scattering  of  the  Og  component  of  the  0  atom  beam.  The  elastic 
scattering  of  the  halogen  molecules  may  be  measured,  but  the  scattering  comes  from 
the  02  as  well  as  0  in  the  0  atom  beam,  and  the  kinematics  for  interpreting 
this  scattering  is  particularly  poor. 

A  crude  estimate  of  the  total  cross  section  may  be  obtained  from  comparing 

the  relative  reactive  and  elastic  signal  intensities  for  these  reactions  with 

18 

the  halogen  atom-halogen  molecule  exchange  reactions.  This  comparison  indicates 

the  cross  sections  are  the  same  order  of  magnitude  for  the  two  reaction  families, 
°2 

roughly  about  10  A  . 

STATISTICAL  COMPLEX  MOOELS 

A  center  of  mass  angular  distribution  which  is  symmetric  about  0  =  90° 
is  suggestive  of  a  reaction  which  proceeds  through  a  collision  complex,  but  is 
not  a  sufficient  condition.  If  a  complex  persists  for  several  rotational  periods 
and  many  vibrational  periods,  all  "memory'1  of  the  initial  conditions  of  complex 
formation  are  lost.  Only  total  energy  and  total  angular  momentum  are  conserved. 
This  situation  permits  simple  statistical  models  to  be  substituted  for  detailed 
studies  of  collision  dynamics.  Many  statistical  approaches  have  been  developed 

Al  l  oo 

including  phase  space  models*^  and  transition  state  models.  ’  These  models 
are  descendants  of  the  transition  state  (RRKM)  treatment  of  unimolecular  decom- 
position  q  and  the  compound  nucleus  treatment  of  nuclear  fission.  Approximate 
transition  state  models  will  be  briefly  discussed  and  applied  to  the  reaction 
systems  under  study  here. 


Approximate  Product  Velocity  Distribution.  The  total  energy  available  to 
the  reaction  products  Is 

Etot  *  E  +  W  +  AD0  »  E1  +  W  (2) 

where  E  and  W  denote  the  initial  relative  translational  energy  and  Internal 

excitation  of  the  reactants,  E'  and  W'  are  the  corresponding  quantities  for 

the  products  and  ADQ  is  the  reaction  exoergicity.  The  role  of  the  statistical 

models  is  to  partition  this  total  energy,  weighted  by  the  density  of  energy 

levels  or  states,  among  the  degrees  of  freedom  of  the  dissociating  complex. 

1 3 

According  to  the  simplest  transition  state  model  ,  the  translational  energy 
distribution  is  given  by 

P(E')  =  A(E')N+(Etot  -  E')  (3) 

The  N+  factor  is  the  energy  level  density  of  "active"  vibrations  and  rotations 
at  the  critical  configuration  of  the  complex.  The  A  factor  takes  into  account 
the  centrifugal  effects  associated  with  complex  formation  and  decomposition. 
The  classical  approximation, 

NXot-  E'>  ■  <Etot-  E'>"  <4) 

is  appropriate  here  because  the  zero  point  energies  are  small  and  similar  for 
reactants,  complex,  and  products.  The  exponent  n  is  given  by  n  *  s  +  hr  -  2, 
where  s  and  r,  respectively,  denote  the  number  of  vibrational  modes  and  active 
rotations  at  the  transition  state.  The  number  of  modes  assigned  to  s  and  r 
depends  on  the  geometry  and  bonding  assumed  for  the  transition  state.  For  a 
"tight"  (all  atoms  bond  in  a  single  "molecule")  three  atom  complex,  a  linear 
structure  has  s  =  4,  r  =  0,  and  n  s  2;  a  nonlinear  one  has  s  =  3,  r  =  1,  and 
n  =  3/2.  In  the  "loose"  Comdex  the  bending  modes  go  over  to  free  rotation  of 


the  product  diatomit  molecule;  thus  s  =  1,  r  =  2,  and  n  =  1. 

In  general  the  calculation  of  the  A  factor  is  quite  complex  because  the 
details  of  the  angular  momentum  coupling  must  be  analysed.  The  magnitude  and 


-  13  - 


direction  of  the  total  angular  momentum  Is  conserved  In  the  collision; 

L  +  J  *  J  *  L'  +  J'  (5) 

where  L  and  L'  are  the  angular  momenta  associated  with  centrifugal  motion  of 

the  reactants  and  products  (L  *  yvb  where  y  is  the  reduced  mass,  v  is  the  velocity 

and  b  is  the  impact  parameter)  and  J  and  J'  the  rotational  momenta  of  the  reactant 

and  product  molecules,  respectively.  To  avoid  consideration  of  the  convolution 
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of  the  unknown  component  momenta  distributions,  the  simplest  approach  is  to 
assume  that  the  centrifugal  angular  momentum  dominates: 

(6) 

As  discussed  later,  this  is  at  least  valid  for  the  entrance  channel.  If  the 
transition  states  for  formation  and  decomposition  of  the  complex  are  taken 
as  the  centrifugal  barriers  in  the  entrance  and  exit  channels,  respectively, 
and  if  these  barriers  lie  sufficiently  far  out  to  be  governed  by  the  r~®  van  der 
Waals  attraction,  then 

A(E')ME7Bm')2/3.  (?) 

“  1  •  E'>V 

The  parameter  B  '  is  the  maximum  exit  centrifugal  barrier  and  can  be  calculated 
if  C  and  C',  the  van  der  Waals  force  constants  in  the  entrance  and  exit  channels 
are  known, 

Bm'  =  (y/y ' )3/2  (C/C* Et.  (8) 

Here  y  and  y'  are  the  reduced  masses  of  the  reactants  and  products.  In  Eq.  (4) 
n  has  been  set  equal  to  2  corresponding  to  tight  linear  complexes.  The  dotted  curve 
in  Figure  12(c)  compares  tfiis  result  for  the  0  +  Br^  reaction  with  the 
experimentally  derived  distribution;  Figure  12(b)  compares  the  corresponding 
product  velocity  distributions. 

Product  Angular  Distribution.  In  the  angular  momentum  limit  given  by 


Eq.  (6),  the  CM  angular  distribution  Is  given  by  the  simple  expression 

I(0)  *  zAme  <9> 

The  curve  given  by  £q.  (9)  is  compared  with  experiment  in  Fig.  10a;  it  is 
clearly  much  too  strongly  peaked  even  with  allowance  for  the  imperfect 
experimental  angular  resolution.  The  less  sharply  peaked  experimental  curve 
must  result  from  the  failure  of  Eq.  (6)  at  least  in  the  exit. channel .  In  the 
entrance  channel  (complex  formation)  the  centrifugal  angular  momentum  must  indeed 
dominate  because  the  rotational  excitation  of  the  reactant  halogen  molecule 
is  very  small  due  to  the  rotational  "cooling"  which  takes  place  in  the  expansion 
of  the  nozzle  beam.  Detailed  consideration  of  the  energetics  and  angular 
momentum  conservation  in  the  product  channel,  however,  indicate  that  the 
rotational  angular  momentum  of  the  diatomic  product  cannot  be  neglected. 

Equation  (6)  must  be  replaced  with 

L%  J  =L'  +  J'  (10) 

As  a  result  the  final  relative  velocity  of  the  products  is  not  constrained  to 
be  perpendicular  to  J.  The  peaking  predicted  by  Eq.  (9)  must  then' be  rounded 


2?  2fi 

A  transition  state  treatment  has  been  developed  *  to  provide  a  simple 
algorithm  for  calculating  the  shape  of  the  angular  distribution.  Three  para¬ 
meters  are  required:  L^,  the  maximum  orbital  angular  momentum  with  which  the 
complex  car  be  formed  (at  the  given  collision  energy);  MQ,  the  root-mean-square 
projection  of  the  total  angular  momentum  on  the  initial  relative  velocity  vector; 
and  M0',  the  corresponding  projection  on  the  final  relative  velocity  vector.  If 

the  ratios  L  /M  '  and  M  /M  '  are  specified,  the  angular  distribution  can  be  read 
mo  oo 
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from  reduced  plots  prepared  by  Miller.  Application  of  this  model  for  the  angular 
distribution  to  these  chemical  systems  is  inconsistent,  however,  since  equation  (10| 


dictates  that  MQ  e  0  and  M  '  t  0.  The  tight  linear  complex  assumed  in 

calculating  the  product  velocity  distribution  also  implies  that  M  '  must 

equal  zero.  Thus,  applying  this  traditional  transition  state  model  would  be 

essentially  an  empirical  fit  of  the  model  to  the  experimental  data.  This  is 

the  approach  that  was  followed  in  the  initial  communication  of  this  work.  The 

0  +  Br^  angular  distribution  was  nicely  fit  by  the  model  when  the  ratio  Ln/M()' 

4 

was  set  at  approximately  2.  Another  inconsistency  is  that  in  order  to  calculate 

l_m  and  the  A  factor  in  the  product  translational  energy  distribution  it  was 

assumed  that  the  formation  and  decomposition  of  the  complex  is  governed  by  the 

centrifugal  barrier  associated  with  the  van  der  Waals  attraction.  This 

means  that  at  the  transition  state  the  complex  must  correspond  to  an  atom  and 

a  diatomic  molecule  interacting  only  through  weak,  non-di recti onal  forces.  Clearly 

this  corresponds  to  a  "loose"  complex  rather.than  the  "tight"  linear  complex 

assumed  in  deriving  the  product  velocity  distribution. 

Loose  Transition  State  Model.  To  avoid  the  problems  discussed  in  the 

preceding- section,  the  simple  models  assuming  a  tight  complex  have  been 

abandoned.  A  model  for  a  three  atom  "loose"  complex  constrained  only  by  Eq.  (10) 

can  be  used  to  obtain  both  the  product  velocity  and  angular  distributions.  This 

model  still  assumes  that  the  critical  configurations  are  the  centrifugal  barriers 

biit  considers  in  detail  the  conservation  of  both  energy  and  angular  momentum. 

24 

The  final  result  is  a  classical  analog  of  the  quantum  phase  space  theory;  every 
state  of  the  complex  at  the  critical  configuration  that  conserves  energy  and 
angular  momentum  is  weighted  equally.  If  the  van  der  Waal's  force  constants  can 
be  obtained  and  if  the  structural  parameters  of  the  reactant  and  product  diatomics 
are  known,  the  model  contains  no  adjustable  parameters.  All  observable  properties 
of  the  reaction  may  be  predicted  including  product  velocity  and  angular  dis¬ 
tributions,  total  reaction  cross  sections,  and  branching  ratios  of  complex  de- 

28 

composition  to  reactant  channel  and  product  channel (s).  This  model  will  be 

16 


described  fully  elsewhere. 
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Table  II  contains  the  parameters  used  In  the  model  calculations.  The 

van  der  Waals  force  constants  for  the  r”®  Interaction,  C,  are  obtained  from  the 

29 

Slater-Kirkwood  approximation.  The  effective  number  of  electrons  were  taken 

to  be  6  for  the  0  atom,  14  for  halogen  molecules,  7  for  halogen  atoms,  and 

13  for  halogen  monoxide  diatomic  molecules.  The  polarizability  values  were 

30 

chosen  to  be  consistent  with  literature  values  and  atomic  additivity.  The 

°3  31 

polarizabilities  (in  A  )  were  taken  as  0.77  for  the  oxygen  atom  and  2.3,  3.1, 

30 

and  5.2  respectively  for  the  Cl,  Br,  and  I  atoms.  The  value  for  each  diatomic 

molecule  was  obtained  by  summing  the  polarizability  of  each  substituent  atom. 

The  induction  terms  were  calculated  using  dipole  moments  (in  D)  of  1.24  for 

IC133,  0.9  for  IBr33,  and  1.26,  1.61,  and  2.45  respectively  for  CIO,  BrO,  and 
34 

10.  The  reactant  internal  excitation  W  was  calculated  from  the  beam  temperatures 
in  Table  I  assuming  classical  rotation  and  that  in  the  nozzle  expansion  no 
vibrational  relaxation  occurs,  but  that  80%  of  the  rotational  energy  is  con¬ 
verted  to  translational  energy.  (This  assumption  is  consistent  with  the  measured 
velocity  distributions  of  the  reactant  beam.)  The  bond  energy  0Q  for  each 

or  or 

molecule  was  taken  from  standard  references  *  with  the  exception  of  10  which 

will  be  discussed  in  a  later  section.  The  diatomic  spectroscopic  vibrational 

constants  uk  and  rotational  constants  B^  were  taken  from  Ref.  36. 
e  e 

As  shown  in  Fig.  10a,  the  angular  distribution  calculated  at  a  nominal 
collision  energy  of  2.11  kcal/mole  gives  excellent  agreement  with  the  CM  angular 
distribution  obtained  by  the  nominal  transformation  of  the  experimental  data. 

The  anisotropy  ratio  B  (equal  to  the  ratio  of  intensity  at  6  =  0°  or  180°  to 
tnat  at  9  =  90°)  predicted  by  the  model  is  2.6  as  compared  to  2.3  obtained  from 
the  nominal  LAB  -*■  CM  transformation.  The  angular  distribution  obtained  in 
another  study^at  a  lower  reactant  translational  energy  (0.8  kcal/mole)  indicated 
a  smaller  value  for  8.  This  is  in  qualitative  accord  with  the  model  which 
predicts  8  *  1.9  at  0.75  kcal/mole. 
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In  contrast  to  the  good  agreement  found  for  the  angular  distribution, 
however.  Figs.  10(b)  and  10(c)  indicate  that  the  predicted  product  velocity 
distribution  Is  in  poor  agreement  with  experiment.  The  loose  transition-state 
model  predicts  too  large  a  fraction  of  the  total  available  energy  entering 
product  translational  energy. 

Comparison  With  Laboratory  Data.  Comparison  of  the  model  predictions  with 
the  curves  derived  from  the  nominal  CM  to  LAB  transformation  may,  however,  be 
misleading.  The  models  are  formulated  in  the  center  of  mass  system  for  a 
complex  of  given  total  energy  formed  in  a  collision  with  specific  relative 
translational  energy.  The  experiments,  however,  were  run  in  crossed  molecular 
beams  and  the  0  atom  beam  was  not  velocity  selected  so  a  wide  range  of  initial 
translational  energies  was  present  in  the  experiment.  Of  less  importance  is 
the  range  of  the  total  energy  of  the  complex  with  respect  to  separated  products 
caused  by  this  spread  in  translational  energies  as  well  as  a  small  range  of 
internal  excitation  in  the  reactant  diatomic  molecule.  Since  no  method  presently 
exists  to  deconvolute  the  angular  and  energy  distributions  from  the  energy 
spreads  and  to  unambiguously  transform  the  data  to  the  CM  system,  a  different 
approach  will  be  taken.  The  model  predictions  will  be  averaged  over  the  initial 
conditions  and  then  transformed  to  the  LAB  reference  frame  to  be  compared  directly 
with  experiment.  Formally,  this  transformation  may  be  expressed  as 

00  oo 

*LAB^6’-)  dV1  S  dv2  W  n2^v2^  (v2/u2)*V*V(V)  ICM(0,u;V)  (11) 
o  o 

Here  (0,v)  is  the  flux  of  the  observed  product  of  laboratory  angle 

and  velocity  x-  *cm  (O.U'.V)  is.  the  center  of  mass  product  flux  predicted  at 

CM  angle  0  and  velocity  u.  This  CM  flux  is  a  function  of  the  relative  velocity 

2  2  h 

of  reactant  approach  V  =  (v^  +  v^  )  when  angular  divergence  of  the  beams  is 

neglected.  The  velocity  number  distributions  in  the  beam  are  given  by  n^(v^) 


and  r^tv^).  The  Jacobian  for  the  CM  to  LAB  transformation  Is  simply  v2/jj2. 

Finally,  V(V)  gives  the  dependence  of  the  total  cross-section  for  complex 
formation  Q  on  the  relative  velocity. 

V 

The  function  I(<^u;V)  may  be  calculated  from  the  loose  transition-state 
model.  Here  for  convenience  an  approximation  is  made.  The  assumption 

I(0,u;V)  *  1(0  ;V)  U  (u;V)  (12) 

is  made  that  the  CM  velocity  distribution  is  independent  of  scattering 
angle  as  given  by  (12).  Except  at  the  limit  of  L  =  J  =  L'  this  cannot  be 
rigorously  true,  but  it  is  expected  to  be  a  good  approximation  for  these  studies. 
This  expectation  is  strengthened  by  inspection  of  Fig.  10(a)  where  this 
separability  appears  to  hold  within  the  uncertainty  of  the  data.  It  has  also 
been  assumed  that  I(G,]j;V)  is  independent  of  internal  state  of  the  reactant 
diatomic  molecule  over  the  small  range  present  in  the  experiment.  Average  values 
for  the  internal  vibrational  and  rotational  energies  can  then  be  assigned. 

The  dependence  of  Q  on  relative  velocity  is  a  direct  consequence  of  the 

V 

loose  transition-state  model.  The  total  cross  section  is  determined  by  the 
largest  impact  parameter  that  still  allows  the  colliding  species  to  surmount 
the  centrifugal  barrier;  thus 

Qc  =  3/2  rr  (2C  /E)1/3  (13) 

tor  the  assumed  van  der  Waals  interaction  which  gives 

V(V)  *  V‘2/3  (14) 

to  substitute  into  Eq.  (11). 

The  hi^oeen  beam  produced  from  the  nozzle  expansion  has  a  very  narrow  velocity 

distribution;  therefore  averaging  over  the  velocity  distribution  of  this  beam  is 
37 

not  nectfsary.  Finally,  no  attempt  is  made  to  calculate  the  absolute  magnitude 

of  the  laboratory  scattering.  The  final  calculation  will  be  normalized  to  the 
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laboratory  data.  With  these  approximations  Eq.  (11)  simplifies  to 

oo 

*LAB^®*— ^  *  5dVl  'W  V^3  ^eiv)  u(uiv)  (15) 

o 

The  calculation  is  accomplished  by  selecting  nine  velocities,  equally  spaced  over 
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the  velocity  range  of  the  oxygen  atom  beam.  At  each  of  these  velocities 
I(0;V)  and  (u;V)  are  calculated.l^(0,v)  is  then  calculated  for  all  desired  product 
laboratory  angles  and  velocities  and  then  summed  over  the  nine  reactant  beam 
velocities.  The  results  for  the  loose  transition-state  model  are  given  for 
the  angular  distributing  by  the  solid  curves  in  Figs.  3 and  4,  and  for  the  velocity 
distributions  by  the  solid  curves  in  Figs.  5-8.  it  should  be  emphasized  again 
that,  except  for  the  single  normalization  constant,  no  adjustable  parameters 
are  included  in  the  calculations. 

Figure  11(a)  compares  the  angular  distributions  predicted  by  the  model  for 
the  five  reaction  systems.  These  distributions  are  for  the  oxygen  atom  velocity 
which  gave  the  greatest  contribution  to  the  total  reactive  scattering  in  the 
integration  over  the  oxygen  atom  velocity.  The  corresponding  collision  energies 
were  3.92,  2.11,  2.44,  2.43,  and  2.41  kcal/mole  for  the  0  +  Cl ^ ,  B^,  and  I^,  IC1 
and  IBr  reactions,  respectively.  The  anisotropy  ratios  calculated  for  the  0  +  I« 
and  IC1  reactions  are  1.9  and  2.8  respectively,  which  are  in  good  agreement  with 
the  respective  empirical  values  of  1.7  and  2.5  obtained  in  Ref.  7.  The  values 
of  Grice  et  al^  are  lower  for  two  reasons:  those  experiments  were  run  at  lower 


collision  energies  where  the  anisotropy  ratio  is  expected  to  be  lower  and  the 
convenient  mathematical  function  used  in  Ref.  7  was  a  two  term  Legendre  polynomial 
expansion  which  rounds  the  CM  #nqular  distribution  at  0  =  0°  and  180°  and  thus 
underestimates  the  anisotropy  ratio. 

The  CM  to  LAB  calculations  were  repeated  using  the  same  angular  distributions 
from  the  loose  transition-state  model  but  using  the  velocity  distributions  corres¬ 
ponding  to  the  energy  distributions  from  Eq.  (3)  for  a  tight  complex  (n  =  2). 


Figure  11(b)  compares  these  energy  distributions  for  the  five  reaction  systems. 

The  results  of  these  calculations  are  shown  by  the  dashed  curves  in  Figs.  3,  4 
and  6-8. 

DISCUSSION:  DYNAMICS 

0  +  Br2  and  0  +  ^  Reactions.  The  laboratory  angular  distributions 

calculated  from  the  statistical  complex  models  are  in  at  least  qualitative 

agreement  with  the  data  for  these  two  reactions.  The  angular  distributions 

appear  in  the  correct  angular  range  and  have  shapes  similar  to  the  experimental 

distributions.  The  calculations  for  the  loose  collision  complex  are  not  peaked 

sharply  enough,  but  the  calculations  for  the  self  contradictory  model  composed 

of  the  angular  distribution  for  a  loose  complex  and  the  velocity  distribution 

for  a  tight  complex  give  nearly  quantitative  agreement.  The  velocity  distributions 

for  the  0  +  Br^  reaction  show  the  same  results:  the  "hybrid"  model  gives  results 

that  agree  in  shape  and  in  position  with  experiment  while  the  loose  collision 
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complex  predicts  too  much  high  velocity  product.  The  agreement  of  the  model 
calculations  with  the  LAB  data,  with  the  nominal  CM  distributions  for  0  + 
and  with  the  anisotropy  ratios  extracted  in  an  independent  investigation^  clearly 
she*  that  the  angular  distributions  are  predicted  very  nearly  exactly  by  the 
loose  collision  complex  model  while  the  velocity  distributions  are  accurately 
given  by  the  approximate  tight,  linear  complex  model.  Grice  and  coworkers^ 
also  found  that  a  tight,  linear  complex  best  fit  their  0  +  I«  velocity  analysis 
data  and  was  in  accord  with  their  0  +  data.  The  angular  distributions 
obtained  by  Grice  et  al^  for  the  0  +  Br^  and  0  +  I~  reactions  are  very  similar 
to  t^ose  obtained  in  this  stud^even  though  their  initial  translational  energy 
was  somewhat  lower  ^  although  our  angular  distribution  for  the  0  +  reaction  is 
better  resolved  than  that  of  Grice  et  al.^  Comparison  of  velocity  spectra  is 
not  possible  as  we  did  not  obtain  this  data  for  the  0  +  I~  reaction  and  Grice 
et  al^  did  not  obtain  velocity  spectra  for  the  0  +  reaction. 
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The  apparent  contradiction  between  the  angular  distribution  and  velocity 
distribution  models  cannot  be  easily  avoided  by  assuming  that  the  product 
velocity  distribution  is  "fixed"  early  in  the  exit  channel  when  the  Incipient 
products  are  still  Interacting  strongly  (tight  complex)  while  the  angular 
distribution  is  determined  at  the  exit  centrifugal  barrier  where  a  “loose" 
complex  must  be  the  proper  description.  In  order  to  calculate  the  angular 
distribution  from  the  loose  complex  all  values  of  L'  (product  orbital  angular 
momentum)  consistent  with  the  conservation  laws  are  weighted  equally.  Since 
L'  Is  related  to  final  relative  product  velocity  V'  by 

L'  =  y'V'b' 

where  b'  is  the  final  asymptotic  impact  parameter  of  the  departing  products, 
only  the  velocity  distribution  from  the  loose  collision  complex  model  is 
consistent  with  that  angular  distribution  unless  some  ad  hoc  assumption  con¬ 
cerning  the  distribution  of  b'  is  made.  Clearly  a  more  unified  model  is 
desirable,  but  for  these  systems  the  angular  and  the  velocity  distributions  can 
be  adequately  calculated  from  the  separate  and  contradictory  models. 

These  experiments  are  particularly  sensitive  to  products  produced  with 
low  velocities  in  the  center  of  mass  system.  This  arises  because  the  Jacobian 
factor  given  in  Eq.  (1)  greatly  magnifies  this  scattering  when  it  is  observed 
in  the  laboratory  system.  The  approximate  transition-state  model  predicts  that 
these  reaction  systems  will  give  a  particularly  large  amount  of  low  velocity 
products  as  may  be  seen  from  inspection  of  Eq.  (7)  and  (8).  The  ratios  y/y' 
and  C/C  are  both  much  smaller  than  unity,  so  the  energy  distributions  peak  at 
very  low  values  as  shown  in  Figure  11(b).  This  effect  is  responsible  for  the 
lack  of  any  clearly  bimodal  LAB  angular  distributions  that  usually  characterize 
reactions  that  proceed  through  collision  complexes.  Even  though  the  C/C'  ratio 
enters  Eq.  (8)  only  to  the  H  power,  the  Jacobian  factor  so  magnifies  the  low 


velocity  scattering  that  the  data  Is  quite  sensitive  to  the  value  of  this  ratio. 
The  fact  that  the  ratios  calculated  fcr  these  experimental  systems  do  give 
adequate  fits  to  the  data  indicates  corroborates  our  choice  of  the  critical 
configurations  for  complex  information  and  decomposition  as  being  the  centrifugal 
barriers  where  the  interaction  between  reactants  or  products  is  adequately 
given  by  the  van  der  Waals  interaction. 

There  are  two  important  factors  which  vary  with  relative  collision  energy 
in  the  model  calculations.  These  are  the  cross  section  for  complex  formation 
and  the  complex  decomposition  branching  ratio  between  reactant  and  product 
channels.  At  low  collision  energies  the  complex  formation  cross  section  is 
larger  (see  Eq.  (13).)  The  fraction  of  complexes  decomposing  into  the  exoergic 
product  channel  rather  than  reforming  reactants  is  also  greater.  Low  energy 
collisions  give  scattering  near  the  halogen  beam,  near  9  =  90°  in  the  LAB. 

If  either  of  these  factors  is  neglected,  the  model  predicts  that  the  angular 
distribution  would  have  nearly  the  same  shape  but  would  appear  at  smaller 
laboratory  angles  than  the  data  show.  These  results  also  give  more  confidence 
in  the  statistical  complex  model  in  general  and  the  centrifugal  barrier  as  the 
critical  configuration  in  particular. 

0  +  Cl 2  Reaction.  For  this  reaction,  the  agreement  between  our  model  and 
experiment  is  not  nearly  as  quantitative.  The  angular  distribution  appears  to 
be  centered  at  snalle*'  laboratory  angles  than  predicted  by  the  statistical  complex 
models.  This  feature  could  be  taken  to  indicate  that  collisions  of  high  velocity 
reactants  which  give  scattering  at  smaller  laborabory  angles  are  more  likely  to 
yield  reaction  than  low  velocity  collisions.  This  would  be  consistent  with  the 

r 

experimental  studyba  which  found  that  this  system  has  an  activation  energy  of 

2.9  kcal/mole.  To  quantitatively  verify  this,  the  calculations  for  the  "hybrid" 

model  were  repeated  with  the  inclusion  of  an  approximate  treatment  of  an  energy 

o 

barrier  in  the  entrance  channel.  At  a  reactant  separation  of  3.4  A  (roughly 
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the  distance  expected  from  the  atomic  radii),  only  the  collisions  with 
translational  energies  along  their  line  of  centers  greater  than  2.9  kcal/mole 
were  allowed  to  form  complexes.  It  was  not  possible  to  easily  account  for 
complexes  decomposing  to  reactants  so  the  variation  of  branching  ratio  with 
collision  energy  was  neglected.  This  approximate  treatment  gave  the  dash-dot 
curves  in  Figs.  (31  and  (6).  The  predicted  angular  distribution  is  in  better 
agreement  with  experiment  insofar  as  the  general  position  is  concerned,  but  it 
is  still  not  peaked  sharply  enough.  The  predicted  velocity  distribution, 

however,  shows  too  much  high  velocity  product  to  agree  with  experiment. 

7d  e 

The  0  +  Cl ^  reaction  has  been  studied  by  Grice  et.  al  ’  at  two  different 
initial  translational  energies,  nominally  3.0  and  6.9  kcal/mole,  using  a  seeded 
oxygen  nozzle  beam.  These  workers  find  that  the  reaction  proceeds  via  an 
osculating  complex^  with  a  ratio  of  backward  to  forward  CM  peaks  of  0.30  at 
the  higher  translational  energy  and  an  increased  ratio  of  0.55  at  the  lower 
energy.  Our  initial  translational  energy  (^2.9  kcal/mole)  is  comparable  to 
their  lower  energy  study  and  coupled  with  the  presence  of  the  modest  activation 
energy  the  angular  distributions  should  be  quite  similar.  Although  they73 
did  not  publish  a  LAB  angular  distribution,  our  LAB  result  would  be  quite 
similar  to  theirs  in  that  we  do  show  a  larger  forward  peak  as  compared  to  the 
backward  peak.  As  we  did  velocity  analysis  at  only  one  angle  (forward  peak) 
we  could  not  construct  a  full  CM  map  and  thus  cannot  compare  directly  with 
Grice  et  als  CM  angular  distribution.  We  do  note  that  they  only  obtained 
velocity  analysis  at  two  angles  at  the  lower  energy  which  did  not  give  a  full 
CM  map  either.  Grice  et  al7e  found  that  there  was  a  dependence  of  CM  velocity 

if 

distribution  on  scattering  angle  and  suggested  that  two  components  were  present, 
one  proceeding  by  a  stripping-like  process  and  one  proceeding  via  the  osculating 
complex.  Our  study,  in  agreement  with  the  work  of  Grice  et  al7c*  at  the  higher 


collision  energy,  predicts  too  much  high  velocity  products.  This  Is  somewhat 
surprising  for  a  system  proceeding  by  an  indirect  process  with  a  more  direct  proc 
superimposed  on  it.  One  would  expect  to  see  more  high  velocity  product  whereas 
we  find  the  opposite.  Although  it  is  likely  that  the  reaction  proceeds  via 
an  osculating  complex,  the  quantitative  determination  of  the  details  of  the 
0  +  Cl g  reactive  scattering  in  the  threshold  region  awaits  further  experiments. 

0  +  IC1  and  I  +  IBr  Reactions.  These  two  reactions  are  complicated  by 
the  presence  of  two  potential  product  channels:  10  +  Cl  or  Br  and  CIO  or 
BrO  +  I.  Despite  careful  searches,  only  the  10  product  has  been  observed  in 
each  reaction.  Figure  5  indicates  that  at  least  in  the  case  of  0  +  IC1,  the 
CIO  product  scattering  is  at  least  an  order  of  magnitude  less  intense  relative 
to  the  10  product  than  is  predicted  by  the  "hybrid"  statistical  complex  model. 
Figure  12  shows  the  branching  ratio  as  a  function  of  initial  reactant  trans¬ 
lational  energy  as  predicted  by  the  statistical  complex  model.  (These  results 
are  from  the  loose  complex  model  and  are  taken  as  the  same  in  the  "hybrid"  model.] 
When  only  the  10  product  is  allowed  the  exoergic  product  channel  is  greatly 
favored  ever  dissociation  of  the  complex  to  reform  reactants  at  low  collision 
energies.  At  high  energies  the  ratios  become  more  nearly  equal.  As  discussed 
above,  this  behavior  favors  scattering  at  angles  near  90°.  In  contrast,  if  both 
product  channels  are  allowed  the  fraction  of  complexes  dissociating  to  the  10 
product  remains  almost  constant  with  reactant  translational  energy.  For  the 
0  +  !C1  "ystem,  only  the  10  channel  can  be  allowed  in  order  to  come  close  to 
matching  the  observed  angular  distribution.  Figure  4  indicates  that  the 
exoeriment?!  data  still  lies  at  an  angular  range  too  large  to  be  fit  by  the  model. 
Cither  seme  factor  that  has  noVbeen  taken  into  account  by  the  model  highly 
favors  reaction  at  very  low  collision  energies,  or  some  experimental  problem 
such  as  the  presence  of  van  der  Waals  dimers  must  account  for  this. 
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The  behavior  shown  In  Fig.  12  may  be  qualitatively  understood  In  terms  of 

the  number  of  states  available  to  the  dissociating  complex  in  each  channel. 

At  low  reactant  translational  energies  the  energy  available  to  the  dissociating 

complex  is  dominated  by  the  reaction  exoergicity.  As  more  energy  is  available 

to  the  complex,  more  states  become  available.  Thus,  at  low  collision  energies, 

the  more  exoergic  channel  dominates.  As  the  reactant  energy  increases,  the 

energy  available  for  each  channel  becomes  relatively  more  nearly  equal.  Here 

the  density  of  energy  levels  comes  to  dominate.  Each  channel  is  simply  an  atom 

and  a  diatomic  at  the  critical  configuration.  The  densities  of  rotational  and 

vibrational  states  increase  as  the  rotational  constant  B„  and  the  vibrational 

e 

constant  u>e  decrease.  The  channels  giving  IC1  (reactants), 

10  and  CIO  should  be  favored  in  that  order.  In  summary,  at  low  reactant 
translational  energy,  the  more  exoergic  channel  dominates  but  as  the  energy  increases 
the  channel  giving  the  diatomic  molecule  with  the  largest  density  of  energy 
states  begins  to  dominate. 

Both,  the  lack  of  any  observable  CIO  or  BrO  product  and  the  position  of  the 
10  product  angular  position  indicates  that  only  the  10  channel  is  allowed  for 
these  two  reactant  systems.  The  most  reasonable  explanation  for  this  behavior 
is  the  geometry  of  the  complex  when  all  thr«  .  atoms  are  interacting  strongly. 

This  will  be  discussed  fully  in  the  next  section. 

O  7 

Both  Lee  et  al  and  Grice  et  al  have  published  results  for  0  +  IC1 

scattering  at  low  collision  energies  and  with  velocity  analysis.  The  results 

7  0 

of  Grice  et  al '  are  very  similar  to  our  results  while  the  results  of  Lee  et  al 

who  employed  an  oxygen  atom  nozzle  beam  show  better  resolution  in  the  center- 

of-mass  and  show  a  double  peak  in  the  LAB  distribution.  Our  "extra"  intensity 

near  the  I  Cl  beam  may  be  due  tp  the  presence  of  IC1  dimers  which  could  yield 

such  scattering. ^  The  velocity  results  obtained  by  the  other  groups  have  been 

fit  using  the  RRKM-AM  5  model.  Grice  et  al  tentatively  interpreted  their  data 
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In  terms  of  a  loose  complex  model  rather  than  the  tight  complex  model  we  have 
8 

used.  Lee  et  al  do  not  present  details  of  their  calculations  so  we  cannot 

discuss  this  point  in  relation  to  their  work.  As  Indicated  in  Figure  4  here, 

the  angular  distribution  is  better  fit  by  the  loose  complex  calculation  at 

the  collision  energy  that  pertains  to  the  experiments  reported  here. 

10  Bond  Energy.  Two  divergent  values  may  be  found  in  the  literature  for 

43  44 

the  dissociation  energy  of  the  10  molecule:  M4  kcal/mole  and  n-57  kcal/mole. 

36 

In  this  study  the  reaction  of  0  with  IC1  (DQ  =  49.6  kcal/mole  )  to  give  the 

product  10  has  been  found  to  proceed  at  low  collision  energies.  The  reaction 

of  0  +  CFjI  (Do  for  CFj-I  =  54  kcal/mole^)  to  yield  10  +  CF^  is  also  spontaneous 

46  47 

at  thermal  energies.  ’  For  these  two  reactions  to  proceed  readily  they  must 
be  either  exoergic  or  only  slightly  endoergic,  so  the  reactant  bonds  broken  in 
these  two  reactions  may  be  taken  as  approximate  lower  limits  for  the  10  dis¬ 
sociation  energy.  These  considerations  certainly  preclude  the  44  kcal/mole  value 
and  suggest  that  the  larger  value  is  approximately  the  minimum  value  consistent 
with  the  results  of  the  0  +  CF^I  reaction. 

The  data  obtained  in  this  study  for  the  0  +  IC1  reaction  provide  a  sensitive 
test  for  the  true  value  of  the  10  dissociation  energy.  This  dissociation  energy 
is  one  parameter  of  the  statistical  complex  models.  A  value  of  57  kcal/mole 
gave  the  best  overall  fit  to  the  angular  and  velocity  distributions.  Variations 
of  2  kcal/mole  gave  significantly  poorer  fits.  For  the  0  +  I?  and  0  +  IBr 
reactions  the  same  Do  value  was  used  in  the  calculations  and  gave  good  results, 
but  since  these  reactions  are  more  exoergic,  the  fit  is  not  nearly  as  sensitive 
to  the  value  of  this  parameter.  We,  thus,  suggest  a  value  of  57  +  2  kcal/mole  for 
the  dissociation  energy  DQ  of  fo.  Other  molecular  beam  studies  of  this  reaction^’® 
determined  a  slightly  lower  value:  53  +  3  kcal/mole.  The  slight  disagreement 
between  the  two  molecular  beam  studies  is  perhaps  due  to  the  assumption  of  a  loose 
collision  complex  model  for  the  velocity  distribution  calculations. 
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Total  Reaction  Cross  Sections.  Table  III  gives  Q  ,  the  total  cross 

section  for  complex  formation  predicted  by  the  statistical  complex  model,  and 

Qr,  the  total  cross  section  for  reaction,  which  is  obtained  from  Qc  by 

multiplying  by  the  branching  ratio  predicted  from  the  loose  collision  complex 

o  2 

model.  The  rough  experimental  values  of  QR  obtained  in  this  study  are  MOA. 

The  pre-exponential  factor  of  the  rate  constant  from  the  flow  tube  study  of 

5a  0  2 

the  0  +  Cl 2  reaction  corresponds  to  QR^1A.  These  experimental  values  are 

significantly  lower  than  the  model  predictions.  Such  has  been  the  case  in 
48 

other  studies  of  reactions  proceeding  through  long-lived  complexes  where  it 
has  been  possible  to  estimate  the  magnitude  of  the  total  cross  section  or 
compare  reactive  and  non-reactive  scattering.  This  behavior  is  usually  interpreted 
as  indicating  that  a  large  fraction  of  the  collisions  which  surmount  the 
centrifugal  barrier  still  cannot  form  a  complex  which  is  capable  of  dissociating 
to  form  products.  Grice^  has  suggested  that  for  these  reactions  only  the  complexes 
which  form  with  a  nearly  coll inear  0-X-Y  orientation  are  capable  of  dissociating 
as  products.  As  will  be  discussed  in  the  next  section,  this  is  in  accord  with 
the  most  stable  geometry  of  the  collision  complex  believed  formed  in  these  reactions 

DISCUSSION:  ELECTRONIC  STRUCTURE 

The  experimental  results  demonstrate  that  the  reaction  of  oxygen  atoms  with 
diatomic  halogens  (excluding  F2)  proceed  via  a  long-lived  complex  or  an  osculating 
complex  in  the  case  of  the  Cl ^  reaction.  Such  long-lived  three  atom  complexes 
probably  have  well-depths  at  least  of  the  order  of  20  kcal/mole  based  on  estimates 
from  RRKM  theory.  The  molecules  XOX  are  thought  to  have  a  * ground  state 
with  an  obtuse  bond  angle^  b^sed  on  the  structure  of  C120^  at  first  glance, 
these  molecules  suggest  themselves  as  viable  intermediates.  However,  one  must 
consider  the  spin  and  spatial  symmetry  constraints  placed  on  the  intermediate 
state  by  the  reactant  and  product  states  assuming  that  the  reaction  occurs  on 
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a  single  electronic  surface.  The  correlations  for  the  reaction  O+Xg-^OX+X 

3 

are  summarized  In  Table  IV.  The  triplet  ground  state  of  the  oxygen  atom,  P  , 
requires  the  reactants  to  have  an  overall  triplet  state  and  this  necessitates 
the  formation  of  a  triplet  intermediate  unless  a  spin  transition  occurs.  This 
spin  correlation  eliminates  the  ground  sate  of  X^O,  'Aj,  from  being  the 
complex  sampled  by  this  reaction.  For  three  atoms,  AB2,  we  need 
consider  only  three  approaches,  C^y,  and  All  of  the  complexes  have  at 

least  one  correlation  within  the  triplet  manifold  with  the  reactant  states. 

The  product  states,  due  to  the  presence  of  two  doublets,  correlate  asymptotically 
with  both  singlet  and  triplet  states  of  the  intermediate. 

In  order  to  provide  further  information  about  the  nature  of  the  electronic 
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structure  complex,  orbital  correlations  must  also  be  examined.  Such  correlations 
are  most  easily  seen  for  the  C  approach  of  an  0  atom  to  the  X0  molecule.  This 

00  V  C 

53 

is  shown  by  the  "weeny-orbital"  diagram  in  Figure  13.  We  have  ignored  the 

valence  s  orbital  as  it  is  low-lying  and  will  not  contribute  to  the  correlations. 

3  -  2 

It  is  obvious  that  the  r  state  of  the  intermediate  arising  from  the  n  configu¬ 
ration  will  be  rc-pulsive  as  two  "sigma"  electrons  on  the  atom  are  approaching  a  two 
electron  sioma  bond.  Thus,  for  the  coll  inear  approach  the  states  of  the 
reactants  must  correlate  with  the  triplet  state  of  the  complex.  However,  for 
the  intermediate  to  show  the  proper  correlations,  the  orbital  occupations  of 
the  complex  must  be  •••  (n*)  (0*)  or  •••  (a*)  (tt*)  in  order  to  give  a  n  state. 

This  is  not  the  (n*)^  configuration  that  we  initially  predicted4  and, 

following  us,  that  Grice  predicted,^  would  be  the  orbital  configuration  of  the 
complex  based  simply  on  state  diagrams. Rather  this  latter  configuration  gives  rise 

•t 


to  states  that  correlate  with  the  repulsive  interaction  described  above.  For 
the  C_  approach,  a  similar  correlation  can  be  expected.  The  orbital  correlation 
diagram  for  the  C^v  approach  is  somewhat  more  difficult  to  construct  but  the 
correlations  ran  be  related  to  excited  states  of  the  molecule  X^O.  As  shown  in 
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Figure  13,  three  Initial  orientations  of  the  0  are  possible.  These  orientations 
correlate  with  the  following  triplet  states  of  XOX  which  we  represent  by  their 
electron  promotions  ir*-»  o*(2)  and  o-a*.  The  promoted  n  electron  arise  from 
the  oxygen  not  the  halogens.  We  note  that  one  of  the  states  of  the  ir+  a*  type 
places  two  electrons  in  an  0  atom  orbital  directed  at  the  o-bond  of  X^.  As  in 
the  similar  case  in  the  approach  we  expect  this  state  to  be  quite  repulsive 
leading  to  a  high  barrier. 
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Investigation  of  the  experimental  evidence  on  the  spectra  of  0C12  and 

the  beam  results  led  to  the  elimination  of  a  C2y  geometry  describing  the 

intermediate.  The  absorption  spectra  for  C120  shows  continuous  absorption 

suggesting  that  all  excited  states  lead  to  dissociation.  The  spectrum  shows 

a  sharp  cutoff  at  15,000  cm~^  (43  kcal/mole)  corresponding  to  the  0+Cl2 

dissociation  limit.  No  spectral  features  at  lower  energies  were  observed. 

Although  this  evidence  for  the  absence  of  bound  excited  states  is  not  completely 

conclusive,  it  is  suggestive  of  the  fact  that  the  accessible  C2y  geometries  are 

probably-  too  high  in  energy  to  be  considered  as  a  viable  intermediate.  Further 

evidence  against  the  C2y  geometry  come  from  our  scattering  experiments. For  the 

reaction  of  0  with  both  IC1  and  IBr  only  the  product  10  was  observed  even 

though  formation  of  the  other  product,  0C1  or  OBr  respectively,  would  be  more 

exothermic.  Since  the  reactions  proceed  via  a  long-lived  complex,  the 

products  should  be  distributed  statistically  and  the  more  exoergic  products 

would  be  produced  in  greater  yield  if  they  are  accessible  geometrically.  From 

the  C2v  geometry,  dissociation  to  either  product  would  be  geometrically  feasible 

and  if  the  C?  geometry  was  being  sampled  both  products  would  be  observed  with 
*-v  # 

the  01  being  the  minor  product.  However,  since  we  do  not  observe  the  other 
products  we  can  eliminate  the  C2y  geometry  from  consideration.  These  results 
lead  us  to  the  conclusion  that  the  geometry  of  the  intermediate  is  an  unsymmetric 
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triplet  of  the  form  OXY  with  either  a  linear  or  bent  configuration  (we  cannot 

distinguish  between  these  geometries  from  our  results)  and  an  orbital  con- 

13  3  1 

figuration  using  the  C  notation  of  either  •••  (o*)  (n*  )  or  •••(**)  (a*)  .Further 

we  note  for  the  case  where  the  halogen  Is  unsymmetrlc ,  that  the  Intermediate 

has  the  halogen  with  the  lowest  electronegativity  In  the  central  position  of 

the  complex.  This  follows  directly  from  our  observation  that  only  01  is 

produced  in  the  0+IC1  and  0+IBr  reactions. 

Our  structure  for  the  intermediate  shows  an  Interesting  feature  on  the 
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Walsh  orbital  diagram  for  these  species.  Walsh  predicted  that  the  upper 

orbitals  in  the  general  diagram  are  ordered  •••Tr,7r?oto*  by  Increasing  energy 

while  Mulliken  predicted^  an  ordering  •••n ,o?tt*o*.  For  our  predicted  triplet 

state  with  orbital  occupancy  or  o*1"*3  the  middle  a*  and  it*  orbitals  must  be 

2 

essentially  degenerate.  If  they  are  not  degenerate  then  a  tt  configuration 

*  *  4  * 

arises  if  o  is  lower  than  tt  while  a  n  configuration  arises  for  ir  lower  than 

★ 

o  .  However,  neither  of  these  configurations  is  allowed  as  the  former  arises 

4 

from  a  repulsive  interaction  as  discussed  previously  while  a  tt  configuration 

*  * 

leads  to  a  singlet  state  which  is  not  allowed.  Since  the  a  and  it  orbitals 

★ 

are  predicted  to  have  a  strong  angular  dependence  with  a  increasing  rapidly 

with  decreasing  angle  then  it  is  unlikely  that  the  molecule  is  strongly  bent 

*  * 

since  this  would  break  the  accidental  o  -n  degeneracy  and  lead  to  incorrect 
states  of  the  intermediate. 

Evidence  from  other  studies  supports  our  conclusions  as  to  the  nature  of 
the  complex.  Rochkind  and  Pimentel^  have  identified  an  isomer  of  Cl^O  in  a 
matrix  isolation  study  of  the  photolysis  of  Cl^O.  Using  infrared  methods,  they 
characterized  a  species  Cl* • -  Cl 0  with  the  CIO  being  similar  to  the  diatomic 
and  the  remaining  Cl  being  weakly  bound  to  the  chlorine  of  the  chatomic.  They 
also  suggest  that  the  molecule  is  not  linear  and  that  the  bonding  between 
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chlorines  follows  the  model  for  the  (p-n*)o  bond  developed  by  Spratley  and 
58 

Pimentel.  A  similar  matrix  Isolation  study  of  FC10  has  been  carried  out 

59 

by  Andrews  and  coworkers.  They  find  a  similar  structure  F**«C10  with  the 

F  weakly  bound  to  the  chlorine  end  of  diatomic  CIO  with  this  moiety  being  very 

similar  to  free  CIO.  The  beam  results  suggested  an  osculating  complex  for  the 

0+CI2  reaction  and  we  can  estimate  the  stability  of  the  0C1  Cl  complex  as 

being  ^10  kcal/mole.  In  agreement  with  these  studies  an  asymmetric  form  was 

predicted  from  the  scattering  experiments.  The  second  study  is  in  agreement 

with  our  prediction  that  the  least  electronegative  atom  (Cl  in  this  case) 

occupies  the  central  atomic  position.  For  the  system  0C1  Cl  it  is  apparent 

that  the  second  chlorine  polarizes  CIO  to  form  a  weak  bond  between  the  chlorines. 

For  complexes  involving  heavier  halogens,  the  well  depths  increase  significantly 

as  scattering  typical  of  a  long-lived  complex  is  found.  Further  evidence  for 

the  electronegativity  ordering  of  the  atoms  in  the  complex  comes  from  studies 

on  the  isoelectronic  XY+  trihalogen  cation  where  the  least  electronegative  atom  is 

+  60 

always  found  to  occupy  the  central  site  even  in  the  case  of  CIgF.  These  studies 

always  show  a  bent  cation. Similar  electronegativity  arguments  are  also  found  to  hold 

9  11 

for  other  reactions  of  radicals  with  halogens.’  We  note  that  the  0+F^  reaction 

has  a  significant  barrier.®  In  this  case,  the  least  electronegative  element, 

0,  cannot  have  easy  access  to  the  central  position  (C£V  approach)  and  thus 

a  barrier  is  expected.  It  is  also  likely  that  the  reactions  proceeds  by  a 

direct  mechanism  rather  than  by  a  complex  mechanism.  Using  the  above  results, 

we  can  summarize  our  predicted  structure  for  the  complex  as  being  an  asymmetric 

1  3 

triplet  OXXwith  an  electron  configuration  of  either  ••••(a*)  (it*)  or 
3  1  f 

•  ••  •  (tt*)j(o+)  and  most  likely  the  complex  is  linear  or  slightly  bent. 

We  now  compare  these  complexes  with  other  known  species  in  order  to  place 
them  in  a  chemical  perspective.  If  we  replace  the  0  atom,  by  a  halogen  the 
complex  has  the  structure  X^  and  is  expected  to  be  quite  weakly  bound  as  is 
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found  from  molecular  beam  studies"  of  the  reaction  X+Yg  which  showed  the 

presence  of  osculating  complexes.  Substitution  of  a  halogen  In  our  complex 

by  an  oxygen  yields  XO2  which  for  CIO2  has  two  stable  isomers  ^  0C10 

and  00C1  In  doublet  states.  Indeed  the  latter  asymmetric  isomer  Is  predicted 

62 

to  be  the  more  stable  form.  Substitution  of  another  oxygen  leads  to  ozone 

which  Is  a  ground  state  singlet  with  a  weak  0.-0  bond.  The  trend  in  molecular 

* 

stabilities  for  these  species  follows  the  addition  of  extra  tt  electrons 

beginning  with  ozone.  We  note  that  only  when  there  is  a  full  complement  of 

v*  electrons  as  in  does  the  system  become  quite  weakly  bound.  If  the 

g 

species  is  HX2  rather  than  then  there  is  no  complex  formed  even  though 
the  system  does  follow  the  electronegativity  ordering  rule  that  the  least 
electronegative  atom  occupies  the  central  site. 
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APPENDIX:  MULTI-ANGLE  TIME  NORMALIZATION 
FOR  ANGULAR  DISTRIBUTION  MEASUREMENTS 

In  many  molecular  beam  experiments  the  scattering  Intensity  at  all  angles 

exhibits  a  systematic  variation  with  time.  Usually  one  reference  angle  with 

a  relatively  noise-free  signal  Is  chosen,  and  Its  magnitude  is  periodically 

measured  to  follow  this  time  variation.  In  this  study,  all  the  data  points 

are  relatively  noisy  and  the  measurement  of  each  is  repeated  several  times. 

Therefore,  a  more  sophisticated  time  normalization  treatment  is  required.  The 

procedure  adopted  here  is  (1)  to  calculate  the  most  probable  value  of  the 

mean  of  all  the  raw  data  points  taken  at  each  angle,  (2)  to  calculate  reduced 

data  points  by  dividing  each  data  point  by  the  calculated  mean  of  all  the  data 

points  taken  at  that  angle,  (3)  to  fit  all  of  these  reduced  data  points  to 

a  Legendre  polynomial  expansion  in  time,  (4)  to  obtain  time  normalized  data 

points  by  dividing  each  original  data  point  by  the  value  of  the  Legendre 

expansion  corresponding  to  the  time  that  the  data  point  was  taken,  and  (5)  to 

repeat  steps  (2)  -  (4)  with  the  mean  recalculated  from  the  time  normalized  data 

points  obtained  in  step  (4).  This  iteration  continues  until  the  sum  of  the 

squares  of  the  deviations  of  the  reduced  data  points  from  the  Legendre  expansion 

converges  to  a  minimum.  It  is  considered  important  here  to  obtain  reliable 

64 

confidence  limits  on  the  final  results;  thus  error  propagation  techniques 
are  used  to  generate  statistically  meaningful  error  bars  on  the  final  averaged, 
time  normalized  signals.  A  paragraph  will  now  be  devoted  to  the  calculations 
required  for  each  of  the  above  steps. 

A  single  measurement,  S..,  of  the  reactive  scattering  at  angle  j  is  the 

in¬ 
difference  of  the  modulated  signal  with  the  beam  flag  open  and  closed.  The 

modulated  signal  is  the  difference  between  the  counts  on  the  "on"  scalar, 

Nft,  and  the  counts  on  the  "off"  scalar,  Ng: 

S1j  ■  (NA  -  Voper,  -  <\  -  Vclosed  (A’  > 
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For  each  of  the  four  scalar  readings  In  Eq.  (Al),  the  distribution  of  the 
magnitudes  for  an  infinite  number  of  measurements  of  S. .  is  expected  to  be 

*  J 

a  Poisson  distribution.  Therefore,  the  standard  deviation  of  the  respective 
distributions  may  be  estimated  from  a  single  reading  by  vff.  The  standard 
deviation  a..,  of  the  distribution  of  S..  values  is  thus  approximated  by 

I  J  I  J 

°1j  -  [("a  +  Vopen  +  ("a  +  Vclose/  (A2) 

Ignoring  systematic  time  varia^.ons  of  the  scattering  intensity  for  the 
moment,  the  most  probable  value  of  the  mean  of  the  distribution  of  S. .  that 

*  J 

may  be  calculated  from  the  finite  number  of  experimental  measurements  is 

•  L  •  1.4.  J  68 

given  by  a  weighted  mean: 

Si  =  ?(S.  ./a.i2)/?n/ai  2)  (A3) 

J  i  ij  ij  i  ij 

The  uncertainty  in  S.  is  given  by  the  estimated  standard  deviation,  o . ,  of  the 

J  J 

distribution  of  many  determinations  of  S.: 

°j  ■  n/fo/ojj2)]*  (M) 

Now  the  reduced  data  points  may  be  calculated: 


(A5) 


The  uncertainty  of  each  reduced  data  point,  o..,  must  also  be  calculated: 

*  J 


I 
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The  time  normalization  function  is  a  weighted  least-squares  fit  of  the 
reduced  data  points  to  a  Legendre  polynomial  expansion  defined  by 

F(t)  =  £  anPn  (cos  t)  (A7) 

_  2 

where  0  <  t  <  180.  Each  of  the  reduced  data  points  is  weighted  by  l/o. .  .  It 

^  3 

is  necessary  to  express  the  time  span,  T,  of  the  experiment  on  the  range  of  t. 

This  coordinate  change  is  accomplished  by  t=cT-d  where  c  =  0.98  [180/(Tmax-Tm^n)] 
and  d  =  c  Tm_.  -  ( . 02 ) ( 1 80 ) .  A  maximum  of  ten  polynomials  is  used,  but  the 
expansion  is  terminated  at  a  lesser  number,  m,  if  two  conditions  are  met.  First 
the  F  distribution  test  of  an  additional  term  as  described  in  Reference  64  is 
applied  after  the  calculation  of  each  additional  term  to  determine  if  that  term 
is  statistically  justified.  If  the  probability  is  ^  5%  or  greater  that  a  set 
of  points  with  random  deviations  from  the  m  term  expansion  would  be  fit  by  the 
m  +  1  term  expansion  as  well  as  the  true  data  points  are  actually  fit,  then  the 
first  condition  is  fulfilled.  Second,  the  x  test  of  Reference  64  is  used  to 
test  the  overall  fit  for  m  terms.  If  the  probability  is  ^  5%  or  greater  that 
the  deviations  of  the  s..'s  from  the  polynominal  expansion  are  as  great  as  would  be 

'  J 

expected  for  the  deviations  of  data  points  with  the  given  uncertainties  from 
the  true  time  normalization  function,  then  the  second  condition  is  fulfilled. 

The  first  approximation,  S..',  to  the  time  normalized  data  points  may  now 

*  J 

be  calculated 

V  ■  VF(t)  (A°) 

The  uncertainty  of  each  S..'  is  given  by  o..'  in 

*  J  *  J 


ij 
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where  °p^)’  the  uncertainty  in  the  time  normalization  function,  is  given  by: 

°F(t)2  '  n  pn2  <cos  ‘>  \  W°> 

where  is  the  uncertainty  in  the  determination  of  a  .  Here  the  covariance 
an  n 

of  S. .  and  F(t)  is  small  and  has  been  neglected. 

•  J 

Now  these  four  steps  are  repeated  with  the  S.  values  for  Eq.  (A5)  calculated 
from  the  time  normalized  data  points  from  Eq.  (A8)  and  uncertainties  from  Eq.  (A9) 
substituted  into  (A3).  Repetition  continues  until  the  sum  of  the  squares  of 
the  deviations  of  the  reduced  data  points  from  F(t)  decreases  by  less  than  5% 
between  two  successive  calculations.  For  all  the  data  reported  here,  two 
repetitions  sufficed,  and  there  was  negligible  difference  between  the  results 
of  the  two  iterations.  Figure  14  shows  the  calculated  F(t)  and  reduced  data 
points  for  the  experimental  run  with  the  most  severe  time  variation. 

Figure  15  gives  the  final  time  normalized  data.  The  weighted  averages  from 
Eq.  (A3)  over  the  final  time  normalized  points  are  considered  to  be  the  best 
values  that  can  be  extracted  from  the  original  data.  Two  sets  of  error  bars 
are  shown;  each  indicates  one  standard  deviation  in  the  determination  of  S.. 

J 

The  left  set  is  calculated  from  Eq.  (A4)  and  the  right  set  is  evaluated  from 
the  scatter  in  the  individual  time  normalized  data  points  at  each  angle.  If 
the  left  set  is  systematically  smaller  than  the  right  set,  it  indicates  that 
the  time  normalization  function  is  trying  to  match  the  inherent  scatter  of  the 
data  and  fewer  terms  should  be  used  in  the  Legendre  expansion. 
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FIGURE  CAPTIONS 

Figure  1  Schematic  of  experimental  configuration.  The  components  shown  are 
surrounded  by  a  large  liquid  nitrogen  cooled  shield  and  enclosed 
In  a  large  vacuum  chamber.  The  reactant  beams  cross  at  an  angle  of 
90°.  Scattered  molecules  can  be  observed  over  the  range  0  =  -15° 
to  119°  in  the  plane  of  the  incident  beams. 

Figure  2  Schematic  diagrams  of  beam  sources,  a)  Microwave  discharge  source 

for  oxygen  atom  beam,  b)  Halogen  molecule  nozzle  beam  source. 

Figure  3  Laboratory  angular  distributions  (number  density)  of  reactively 

scattered  products.  The  circles  indicate  the  weighted  averages  of 
all  data  taken  at  each  angle;  the  error  bars  indicate  one  standard 
deviation  (see  the  Appendix  for  complete  discussion.)  The  0  + 
data  are  an  average  of  two  experimental  runs.  The  scattering  angle 
is  in  the  plane  of  the  parent  beams.  The  direction  of  the  0  atom 
beam  is  at  0°,  the  halogen  beam  at  90°.  The  vertical  arrows  indicate 
the  position  of  the  centroid  velocity  vector  calculated  for  the  com¬ 
bination  of  beam  velocities  of  the  most  probable  collision  energy. Curves 
show  distributions  calculated  from  collision  complex  models  using 
the  loose  (solid  curves)  or  tight  linear  (dashed  curves)  transition 
state  models.  An  approximate  calculation  for  0  +  assuming  an 
activation  energy  of  2.9  kcal/mole  is  also  included  (dot-dash  curve). 

Figure  4  Laboratory  angular  distributions  of  reactively  scattered  products. 


Curves  and  symbols  are  defined  in  Fig.  3. 
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Figure  5  Laboratory  angular  distributions  of  reactlvely  scattered  products 
from  the  0  +  IC1  reaction.  The  data  points  Indicate  the  search 
for  the  CIO  product.  Curves  Indicate  the  expected  CIO  product 
(solid  curve)  compared  to  the  10  product  (dashed  curve)  as  predicted 
by  the  loose  transition-state  model.  Equal  detection  efficiencies 
for  the  two  products  were  assumed  to  normalize  the  curves  and 
data  points. 

Figure  6  Laboratory  velocity  spectrum  (flux  density)  of  reactively  scattered 
CIO  product  at  a  scattering  angle  of  60°  for  the  0  +  C^  reaction. 
The  error  bars  indicate  one  standard  deviation  and  include  the 
uncertainty  in  the  normalization  to  the  angular  data  as  well  as 
that  in  the  velocity  spectra  measurement.  The  vertical  arrow  indi¬ 
cates  the  upper  limits  on  the  velocity  from  energy  conservation  for 
the  nominal  collision  energy.  Curves  show  distributions  calculated 
from  the  statistical  complex  models  convoluted  over  the  velocity 
resolution  of  the  detector.  The  solid  curve  is  for  a  loose  complex; 
the  other  two  are  for  a  tight  linear  complex;  the  dot-dash  curve 
includes  an  approximate  treatment  for  an  activation  energy  of  2.9 
kcal/mole. 

Figure  7  Laboratory  velocity  spectra  of  reactively  scattered  BrO  product  at 
indicated  scattering  angles  for  the  0  +  reaction.  The  curves 
(as  in  Fig.  6)  are  calculated  from  the  statistical  complex  models; 
here  convolution  over  the  velocity  resolution  of  the  detector  was 
not  necessary. 

Figure  8  Laboratory  velocity  spectra  of  reactively  scattered  10  product  at 
indicated  scattering  angles  for  the  0  +  IC1  reaction.  The  curves 
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(as  in  Figure  6)  are  calculated  from  the  statistical  complex  models 
with  no  convolution  over  the  velocity  resolution  of  the  detector. 

The  three  horizontal  bars  Indicate  this  velocity  resolution  (FWHM) 
at  the  respective  velocities. 

Figure  9  Polar  (velocity-angle)  contour  map  of  BrO  reactively  scattered  flux 
intensity  in  the  center  of  mass  system.  The  0  atom  beam  direction 
is  0°  and  that  of  the  halogen  molecule  beam  is  180°.  Ticks  along 
the  radial  lines  mark  velocity  increments  of  100  m/sec. 

Figure  10  Product  distributions  in  the  CM  system  for  the  0  +  Br^  reaction. 

Experimental  curves  (solid  lines)  were  derived  from  Figure  9  . 

The  results  of  the  loose  (dashed)  and  tight-linear  (dot-dash) 

transition-state  models  are  shown.  The  angular  distributions  (a) 

are  each  normalized  so  that  the  total  scattering  given  by 

2tt  J*  1(e)  sin  ede  is  equal  to  unity.  Included  in  (a)  is  the  curve 
o 

(dotted)  expected  if  orbital  angular  momentum  dominates.  The  velocity 
distributions  (b)  and  the  corresponding  relative  translational 
energy  distributions  (c)  are  all  normalized  so  that  their  total 
integrals  equal  unity. 

Figure  11  (a)  CM  angular  distributions  calculated  from  the  loose  transition- 

state  model.  The  curves  have  been  arbitrarily  normalized  to  1.0 
at  6  =  0°.  The  0  range  90°  -  180°  is  a  reflection  about  90°  of 
the  0°  -  90°  range,  (b)  Product  translational  energy  distributions 
calculated  from  the  tight  linear  transiton  state  model.  For  clarity 
the  curves  for  the  0  +  IC1  and  0  +  IBr  reactions  have  been  lowered 
by  one  unit  along  the  ordinate.  The  curves  are  normalized  so  that 
their  integral  is  equal  to  unity. 
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Figure  12  Probability  of  complex  decomposition  to  each  channel  as  a  function 
of  Initial  translational  energy  for  the  0  +  IC1  reaction.  Dashed 
curves  pertain  to  when  the  I  +  CIO  channel  is  forbidden,  solid 
curves  to  when  all  three  channels  are  allowed. 

Figure  13  Orbital  correlation  diagram  for  the  C9  and  C  approaches  of  an 

b  V  ®v 

3 

0  atom  (  P  }  to  a  diatomic  homonuclear  halogen  a)  C$  approach. 

The  state  symmetry  is  shown  on  the  left  while  the  orbital  con¬ 
figurations  of  the  complex  are  shown  on  the  right.  Note  the 
repulsive  interaction  for  the  Z  approach,  b)  C^y  approach.  The 
notation  for  the  excitation  corresponds  to  the  type  of  orbital 
excited  from,  in  this  case  a  a  bonding  orbital  or  a  n  out-of-plane 
orbital,  to  an  anti-bonding  a*.  We  have  not  tried  to  identify 
which  orbital  is  the  o  or  n  orbitals,  but  presumably  they  are  the 
6a^  and  2b^  orbitals  respectively.  The  a*  orbitals  would  be 
expected  to  be  the  7a^  and  5b2  orbitals.  Note  the  repulsive  inter¬ 
action  in  diagram  on  the  right. 

Figure  14  Time  normalization  function  and  reduced  data  points  for  the  0  +  1^ 
reaction  data  presented  in  Figure  3.  Each  reduced  data  point  is 
obtained  by  dividing  an  individual  measurement  of  the  signal  at 
a  given  angle  by  the  mean  of  all  the  measurements  made  at  that 
angle.  The  abscissa  scale  is  obtained  by  linearly  normalizing  the 
time  span  of  the  experiment  to  a  0  to  180  span  for  purposes  of  the 
Legendre  polynomial  expansion.  At  t  «  40,  the  0  beam  discharge 
source  was  retuned  resulting  in  a  sharp  increase  in  product  scatter 

tf 

ing  intensity. 
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Figure  15  Time  normalized  data  points  (open  circles)  and  weighted  averages 
(solid  diamonds)  for  data  of  Fig.  14.  The  error  bars  Indicate 
two-thirds  confidence  limits  (one  standard  deviation)  for  the 
determination  of  the  average.  The  left  bars  are  calculated  by 
the  error  propagation  techniques  discussed  in  the  text*  and  the 
right  bars  are  evaluated  from  the  scatter  in  the  individual  data 
points  shown. 
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TABLE  I. 

PARENT  BEAM  PARAMETERS3 

Reactant 

To 

tn  pn 

ci2 

1200 

683 

400 

Br, 

1000 

473 

870 

l2 

1200 

473 

150 

IC1 

1200 

418 

260 

IBr 

1200 

402 

160 

3Tq  and  denote  the  temperatures  (°K)  of  the  oxygen  atom  beam  and  nozzle 
source  respectively.  Tq  was  determined  by  velocity  analyzing  the  oxygen 
atom  beam  under  identical  conditions  and  fitting  the  results  to  a  Maxwell- 
Boltzmann  velocity  distribution  to  extract  the  temperature.  P^  refers  to 
the  halogen  gas  pressure  (torr)  behind  the  nozzle;  this  pressure  is  determined 
from  the  vapor  pressure  of  the  halogen  at  the  temperature  of  the  supply  reservoir 
except  for  Cl ^  for  which  the  pressure  was  approximately  measured  by  a  gauge. 


* 
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TABLE  II.  PARAMETERS  FOR  TRANSITION- STATE  MODELS* 


Parameter  0  +  Cl  2 

0  +  Br2 

0  +  i2 

0  +  IC1 

0  +  IBr 

C(10"12  erg. ft  6) 

96 

117 

165 

135 

142 

-1  9  0  6 

C* (10  ^erg.A  ) 

173 

258 

541 

281 

351 

309 

375 

W(kcal/mole) 

1.0 

0.7 

0.9 

0.6 

0.6 

D0  XY(kcal/mole) 

57.2 

45.4 

35.6 

49.6 

41.9 

D.  X0(kcal/mole) 

63.3 

55.3 

57 

57 

57 

0 

63.3 

55.3 

ECQLL(kcal/mole)b 

3.92 

2.11 

2.44 

2.43 

2.41 

u»e  XY(cm_1) 

559.7 

323.3 

214.5 

384.3 

268.7 

u  XOfcnf1) 

780 

771 

681.4 

681.4 

681.4 

e 

780 

771 

Be  XY(cm‘1) 

0.2441 

0.0811 

0.0374 

0.1141 

0.0561 

B„  XO(cm_1) 

0.622 

0.4278 

0.3403 

0.3403 

0.3403 

e 

0.622 

0.4278 

aFor  the  quantities  where  two  entries  are  given  for  the  mixed  halogen  reactions, 
the  upper  value  is  for  the  10  product  channel  and  the  lower  value  is  for  the 
alternate  (CIO  or  BrO)  product  channel. 

initial  collision  energy 
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TABLE  III.  CROSS  SECTIONS  CALCULATED  FROM  LOOSE  TRANSITION  STATE  MODEL3 


;• 

L 

Reaction 

o2 

QC(A  ) 

o* 

VA ) 

.. 

w 

0  +  ci2 

42 

35 

Si 

► 

0  +  Br2 

55 

45 

0  +  i2 

59 

50 

0  +  IC1 

55 

42 

0  +  IBr 

56 

43 

s 

The  quantity  Qc  is  calculated  from  Eq.  (13)  of  text  at  the  nominal  collision 
•  energy.  QD  is  obtained  by  multiplying  Q  by  the  branching  ratios  predicted 

front  the  loose  complex  model  averaged  over  reactant  translational  energy.  Only 
I  the  10  channel  was  allowed  in  the  mixed  halogen  reactions. 

4 


4 
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TABLE  IV.  STATE  CORRELATIONS  FOR  0+X2 

REACTIONS 

Approach 

Reactants 

Intermediate 

Products 

Geometry 

(Possible) 

[X(2Pu)+0X(2tt)] 

^ooy 

3_-  .  3 

Z  +  u 

1  +  b 
+  2 

3  i  c 

+  V  TT 

3A+3tt+V+3E- 

+1A+1tt+1E++1z' 

Cs 

2(3A')+3A" 

^A' +1a"+3A'  +3A"d 

3(3A')+3(3A") 

+3( 'A' )+3( 'A" ) 

C2v 

7S- 

3a2+3b1+3b2 

1A  e 

♦VVS  f 

+1A]+1B]+1B2  f 

o  n 

2(3A1)+2(3A2)+3B1+3B2 

+2(1A1)+2(1A2)+1B1+1B2 

These  states  arise  from  the  •••(cr*)2(7r*)2  electron  configuration. 


bThis  state  arises  from  the  •••(tt*)4  electron  configuration. 


cThese  states  arise  from  either  the  •'••(tt*)3(o*)1  or  •••(o*)1  (it*)3  electron 
configurations.  The  singlet  states  of  the  C  intermediate  could  arise  from 

0OV 

the  0(  0^)  excited  state  which  gives  ^"Vtt+^a. 

dThese  states  arise  from  the  •••(a')2,  •••(a")2  or  •••(a* J1 (a")1  configurations. 
0 

This  state  arises  from  the  ground  state  configuration. 

fThese  states  arise  from  the  confgurations  ••• •(bj J1 (a^ )] •••(a1 (a] )]  and 
(a-j estimated  from  Walsh  diagrams. 
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